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William M. Adams, Jr. and Sherwood H. Tiffany 
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Hampton, Virginia 


ABSTRACT 


A control law is developed to suppress symmetric flutter for a mathematical 
model of an aeroelastic research vehicle. An impleraentable control law is attained 
by including modified LQG (Linear Quadratic Gaussian) design techniques, controller 
order reduction, and gain scheduling. An alternate (complementary) design approach 
is illustrated for one flight condition wherein nongradient-based constrained optimi- 
zation techniques are applied to maximize controller robustness. 
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upper and lower limits on constraint variable vector 
vector of inequality constraint violations 

state, control and noise distribution matrices in plant output 
equations F 

matrix of viscous damping force coefficients 
regulator gain 

reduced order controller gain magnitude, degrees/g 

gain magnitude of initial, unscheduled, reduced order controller 
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gain margin 

plant transfer function (y/u R ), FS-off 
acceleration due to gravity 
altitude 

controller transfer function (u FS /y) 

regulator and estimator performance indices 

generalized stiffness matrix 
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Kalman estimator gain 

generalized mass matrix 
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number of counterclockwise encirclements of the -1 point by the 
open loop transfer function as co varies from 0 to ~ 

number of generalized coordinates, controls, gusts, outputs, and 
unsteady aerodynamic lag coefficients 

steady-state convariance of state estimate 

matrix of generalized aerodynamic force coefficients and its 
s-plane approximation, respectively, 

dynamic pressure 

intensities of control, gust and measurement noises 
state and control weighting matrices in regulator design 
closed loop, steady-state covariance matrix 
Laplace variable 

fixed portion of reduced order controller 
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negative of commanded control input from FS control law 
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is process and measurement noise intensity 
matrix 

diagonal weighting matrix in augmented performance function (see 
Eq. (22)) 

gust velocity 
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w =(w ,w ,w ) is a vector of zero mean, white noise errors in 
u ’ g’ m * 

the input, white noise driving the gust filter and white noise in 
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plant state vector estimate 
vector of control rotations 

outboard aileron and commanded outboard aileron rotations, 
respectively 

vector of generalized coordinates 
minimum singular value 
phase margin 

augmented performance function minimized in constrained optimiza- 
tion solution for robust controller (see Eq. (22)) 

frequency, rad/ sec 

reduced order 

expected value 



INTRODUCTION 


This paper describes the design of a control law for suppression of symmetric 
flutter for a mathematical model of the DAST (Drones for Aerodynamic and Structural 
Testing) ARW-2 (Aeroelastic Research Wing Number 2) aircraft. An implementable con- 
trol law is attained by including the following elements in the design process: 
development of linear, time-invariant state space evaluation and design models that 
include unsteady aerodynamic force effects (Refs. 1-4); use of LQG (Linear Quadratic 
Gaussian) methods for full order controller design including robustness recovery 
(Refs. 5,6); definition of a candidate reduced order controller (Refs. 7-9); and 
development of a gain schedule to improve off design performance. 

Results are presented which indicate how the evaluation and design models of 
the plant are defined. The control problem is exhibited by displaying the flutter 
characteristics of the uncontrolled vehicle and by stipulating the design criteria. 
A description is given of the methodology employed in determining the full and 
reduced order controllers; and the performance of the scheduled reduced order 
controller is shown both at the design point and at off design conditions. 
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The DAST ARW-2 will employ what is effectively a single-input/single-output 
control law for suppression of flutter in the symmetric degrees of freedom. Conse- 
quently, robustness can be characterized in terms of gain margins, phase margins, 
minimum eigenvalue (singular value) of the return difference transfer function and 
high frequency characteristics of the controller. Robustness constraints were not 
explicitly included in the design algorithm used to define the scheduled reduced 
order controller. Instead, robustness characteristics were examined offline and 
improved by repeating the design using the robustness recovery procedure (Ref. 6). 

One example is also shown where constrained optimization techniques 
(Refs. 10-13) are utilized to explicitly include robustness criteria within the 
design algorithm. The minimum singular value of the return difference matrix is 

maximized subject to constraints on stability, control power, gain margins, and phase 
margins. A similar approach, which employs a gradient-based optimizer, has been 
applied to improve the robustness of a lateral stability augmentation control law 
(Ref. 14). The nonlinear programming algorithm (Refs. 15,16), selected here to 
perform the optimization, requires no gradient computation and has recently been 
successfully applied to several design problems (Refs. 13,17,18). 


DAST ARW-2 SUMMARY 


The primary objective of the DAST ARW-2 research is to provide a partial assess- 
ment of the validity of applying current analysis and design techniques in the design 
of actively controlled aircraft. In this activity a remotely piloted drone which 
will Incorporate several active control functions has been designed and built and 
will be flight tested. Correlation of experimentally measured and analytically 
predicted performance will be documented. 

The ARW-2 wing has a supercritical airfoil, a 10.3 aspect ratio and a 25° sweep 
at the quarter chord. It was purposely designed to require flutter suppression (FS) , 
maneuver load alleviations (MLA) and gust load alleviation (GLA) in some flight 
regimes. Furthermore, the wing was positioned on the fuselage so as to require 
relaxed static stability (RSS) augmentation to achieve satisfactory handling charac- 
teristics. 

An early phase of research (Refs. 19,20) defined actuator requirements, identi- 
fied desirable sensor locations and developed preliminary active control laws. 
Refinement of the modeling and control law design has resulted in the analog imple- 
mentation of the control laws into hardware. This hardware has been constructed such 
that gain and filter constants can be changed and new compensation can be added, if 
necessary, as the knowledge of the aircraft becomes more precise. The undocumented 
structural model employed in this paper is the same as that used in obtaining the 
currently implemented control laws. It is similar to that of references 19 and 20. 

The wings have been constructed at NASA Langley, the control hardware and actua- 
tors have been installed, and ground vibration and static loads tests with the wing 
cantilevered have been conducted. A more detailed structural model is under develop- 
ment which is being tempered by the experimental data. When the updated structural 
model is completed, the FS design will be repeated for the symmetric case and an 
antisymmetric FS law will be developed. 
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Figure 1 depicts the sensors and control surface that will be employed in the 
symmetric FS control law and defines how the sensor signals are separated into 
symmetric components prior to compensation. 


MATHEMATICAL MODELS 


Evaluation Model 

Modal Characteristics. - A modal characterization of the aircraft is employed 
which resulted from performing a free-free vibration analysis of the ARW-2. Twelve 
modes are retained for the symmetric degrees of freedom which are mean-body vertical 
displacement and pitch and ten symmetric elastic modes. Table 1 lists the natural, 
in vacuum, frequencies associated with each of the elastic modes and figure 2 shows 
the normal components of the first six symmetric elastic mode shapes. Fore and aft 
motion, not shown in the plots, predominates for elastic mode 3. 

Unsteady Aerodynamics Force Computation. - Unsteady aerodynamic force coeffi- 
cients are computed using a doublet lattice (Refs. 21,22) module contained in the 
ISAC (Interaction of Structures, Aerodynamics and Controls) program (Ref. 1). The 
paneling of the lifting surfaces that was employed is shown in figure 3. The circles 
shown indicate points at which modeshape data are defined. Aerodynamic force coeffi- 
cients are computed for a number of Mach numbers and, at each Mach number, for the 
following set of reduced frequencies. 

(0.0, 0.005, 0.01, 0.05, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.8, 1.0) 

Reduced frequency satisfies the following relationship 


k 


wb 

U 


(i) 


where b Is a reference length (chosen here to be one half the mean aerodynamic 
chord of the wing), U is airspeed and w is frequency of oscillation in rad/sec. 

Frequency Domain Equations. - Given the modal characteristics and associated 
aerodynamic force coefficients, a frequency domain form of the equations of motion 
can be written. 


[— M(i) 2 +Diu+K+qQ^ ]{ 5 } + [-M 6 <o 2 +qQ fi ^ Q g ]{ * } = {o} (2) 

where 

{?} is a vector of generalized coordinates 

{5} is a vector of control rotations 

{Wg} is a vector of gust velocity components 

M,D,K are generalized masses, dampings and stiffnesses, respectively 

(damping ratios of 0.005 were assumed for each elastic mode) 

is control mass coupling matrix 
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Q(k,M)=[Q ? :Q 6 :Q g ] 


is a matrix of unsteady aerodynamic force coefficients 
which is complex and frequency and Mach number dependent 


and 

9 is dynamic pressure 

In this form of the equations of motion, the assumption is made that control surface 
actuators are sufficiently powerful that aerodynamic hinge moments and inertial cross 
coupling hinge moments do not affect control position (rotation). 

The following transfer function representation of the actuator with aileron 
attached is employed in both the evaluation model and the design model to relate 
control deflection to commanded control deflection: 


6 a = il 4 _? 6 (3) 

a s + 180 s 2 + 251s + 314 2 a c 

This representation fits the experimentally determined frequency response up to 300 
rad/sec. 

S-plane Aerod ynamic Force Approximation .- An approximation to the unsteady aero- 
dynamic force coefficients is made in order to convert equation (2) into a set of 
first order, time-invariant state equations more amenable to linear systems analysis 
and design techniques. A least squares curve fit is made, for a specific Mach 
number, of the matrix of frequency dependent coefficients using a matrix function of 
the form (Refs. 3,23-26). 
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In the evaluation model, n^ is 4 and the b^ are 0.0939, 0.1878, 0.2817, 
0.3756. Linear constraints are imposed upon the rigid body columns of the A-.’ 
matrices which require that the curve fit match the tabular data and its slope at 

k - 0 (Ref. 27). In addition, the A 2 column corresponding to the gust mode is 
constrained to be zero. 


State Space Equations of Motion. — Under the assumption that the curve fit in 
equation (4) for s = i U/b k is valid for points off the imaginary axis, linear 
time-invariant equations of motion can be written of the form 


x = Ax + Bu+Bw 
u w 


with outputs 


where 


y=Cx + Du+Dw 
u w 


(5) 
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The evaluation model plant state vector is 77 by 1 . 


Design Model 

A reduced order model of the plant was developed in order to lower the cost of 

the design of the control law. Only two lag terms were retained in the s-plane 

approximation of the unsteady aerodynamic forces, equation (4). The values selected 
for the denominator coefficients, (b^ = 0.31 ,b 2 = 0.71) yield a good fit to the 
computed aerodynamic data. 

A subset of the modes in equation (2) was selected for retention. Modes with 

jiatural frequencies below or near the open loop flutter frequency (116 rad/sec for 
M - 0.86) which were observed to have little effect upon the flutter characteristics 
were truncated. Vertical translation pitch, and the second elastic mode, (see 

Fig. 2), were truncated. Modal residualization (Refs. 4,8,13,29) was employed to 

retain the static effects of the four highest frequency modes. Thus, five modes were 
retained in the design model. The resulting state space representation of the design 
model of the plant has 25 state variables (10 vehicle, 10 aerodynamic, 3 actuator, 2 
gust) . 


Figure 4 is a Bode plot comparison of the amplitude of the symmetric signal to 
be sent to the controller per unit commanded control input for the design and evalua- 
tion models. This figure shows that the reduced order design model is in good agree- 
ment with the full order evaluation model (phase angle comparisons, not shown, also 
exhibit excellent agreement at frequencies below 500 rad/sec). 
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DESIGN OF CONTROL LAWS 


Design Criteria 

Flutter Bound aries. - Figure 5 shows three boundaries. The predicted flutter 
boundary (FS-off) is shown as the solid line. Flutter occurs, for the uncontrolled 
aircraft, to the right of this line. The dashed line, which will be referred to as 
the nominal stability margin boundary, is a boundary to the left of which the FS 
control law should provide stability with ±6dB gain margins and ±45° phase margins. 
The third curve, denoted as the desired closed loop flutter boundary, defines a 

boundary to the left of which the FS control should provide stability with no 
constraints on stability margins* 

^ on ^ ro ^ LS P es ^g p and Evaluation. Points *— The control law design point in this 
study is at an altitude of 15,000 feet and a Mach number of 0.86. This point is 
indicated on figure 5. Off design performance of the control law is evaluated at 
selected points which are also shown on figure 5. 

Control Con straints .— The peak deflection and rate capabilities for the out- 
board ailerons are ±15° and ±740°/sec, respectively. The control law is to be 
designed such that saturation does not occur, in the root mean square (rms) sense 

for an input gust spectrum having an rms of 12 feet/sec. The control saturation 

constraints are to be satisfied at all points to the left of the nominal stability 

margin boundary of figure 5. The Dryden spectrum is assumed for the design phase and 
the von Karman spectrum is used for controller evaluation. 

Altitude Root^ Loci.- Figure 6 shows the characteristic roots as a function of 
altitude (FS-off, M = 0.86) for the symmetric evaluation model. Flutter is predicted 
to occur at 18,750 feet (q = 531 lb/f t 2). The design point dynamic pressure is 

16.5 percent higher than the dynamic pressure at the flutter point. 


Full Order Controller Design 

The ISAC (Ref. 1) program was used to generate the reduced order design model of 
the plant. The resulting linear, time— invariant equations are 


x = Ax + Bu+Bw 
u w 

- - - ( 6 ) 

y = Cx + D w 
w 

where (A,B^) is stabilizible and (A^ ,C^) is reconstructible. 

LQG algorithms in the 0RACLS (Ref. 5) program were then employed to obtain full 
order controllers, with respect to the design model. Figure 7 depicts the block 
diagram representation of the closed loop system and defines the inputs and outputs 
that will be examined to evaluate the performance of candidate controllers. 

Regulat or Design. - Steady-state regulator gains were determined which minimized 
a performance index of the form 
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(7) 


J R " "t“ E {/ 0 l (x T (t)R 1 x(t)+Up S (t)R 2 u ps (t))dt} 


where E is the expected value operator. Rj is null and R, is the identity. 
This choice of weighting matrices, plus the constraint that the closed loop system be 
stable, results in a solution which reflects unstable poles about the imaginary axis, 
leaving stable poles fixed. This is also the minimum control effort solution which 
stabilizes the system (Ref. 7). Nevertheless, by reflecting the unstable poles, it 
provides a finite margin of stability. 

Estimator Design. - Steady-state Kalman filter gains were determined which mini- 
mize the performance index 


I t, 

J E " lim T E t / 0 e (^(Odt} 


tj*® 1 


( 8 ) 


for the system defined in ( 6 ) where 

e(t) « x(t) - x(t) ( 9 ) 

and x(t) satisfies (The reference to t is dropped for compactness.), 


x = Ax + B u u + L [y-Cx] 
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(ID 


so that the state excitation and observation noises are correlated. The correlation 
is a result of the modeling of aerodynamic forces due to gust inputs (Eq. (4)). The 
coupling term has a minor impact on the resulting control law designs. The solution 
for the steady-state Kalman filter gain is given by (Ref, 30) 


L = [pc t +v 12 ]v ~ 1 


( 12 ) 


Riccatti P equation Steady " State C0Variance of the estimate, satisfies the algebraic 


[a- v 12 V 2 C]P + P[;-V 12 V 2 1 C] T + P^V" 1 ^* + 
The estimator equation may be rewritten as 


(13) 
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(14) 


x = [A - B u F - LC]x + Ly 


where upg - Fx and -F is the solution to equation (7). 

Robustnes s Recovery .- Full state feedback control laws designed using LQ method- 
ology are known to possess good robustness properties (Refs. 31,32). LQG designs, 
however, may exhibit unsatisfactory robustness characteristics. The introduction of 
fictitious process noise at the control input has been shown to provide robustness 
recovery in LQG designs (Ref. 6). Consequently, several estimator designs were 
performed in which the process noise intensity, Ry, was varied to obtain an accept- 
able tradeoff between stability margins, control law bandwidth, and available control 
power. Results for two of the designs are summarized in figures 8 and 9. The loop 
transfer function Nyquist plots of figure 8 show how the stability margins improve 
with increasing process noise at the actuator input (the full order evaluation model 
of the plant was employed in constructing Fig. 8). The Bode plots of figure 9, how- 
ever * reveal the undesirable feature that the commanded control input at high 
frequencies becomes larger with increasing input process noise. The controller 
corresponding to R u = (0.79 deg) 2 was selected as a basis for initiating the 
search for a reduced order controller. Note that this controller meets the specified 
stability margin criteria. Its control requirements, not shown, are also well within 
the ±15°, ±740°/sec capabilities of the actuator. 


Reduced Order Controller Selection 


The 25th order controller developed using LQG techniques was reduced to obtain 
an implementable control law. Figure 10 shows the closed loop block diagram of plant 
and controller where the controller state space representation is given by 


and 


x = A x + B y 
c c c c 


(15) 


u 


FS 


C x 
c c 


(16) 


For the case of a full order controller, A c and B„ are defined by equa- 
tions (13) and (14) and C c = F. C M 


Order reduction was initiated by transforming (15) to block diagonal form and 
examining the poles, zeros and residues of the full order (upg/y) transfer func- 
tion. Modal truncation was then performed in the transformed domain to obtain the 
following candidate reduced order controller. 


U FS F 0 692 * 3 (s+98.59)(s+213.5) s 2 +214.9s+164.5 2 s 2 +l .030s+136.4 2 

y s+692.3 s 2 +118.8s+66.29 2 s 2 +58.48s+94 .42 2 s 2 +2.411s+136.9 2 

b 2 +251 . ls+314 . I 2 
s 2 +406 . 1 s+567 .9 2 
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where F 0 - 0.5064 deg/g was selected so that the static gain of the full order 
and reduced order controllers were the same. 

The sharp notch at 136 rad/sec essentially cancels the peak in the (y/u) trans- 
fer function due to the lowly damped fore and aft wing mode. In some cases, e.g. if 
the flutter mode frequency were nearly coincident with the notch, failure of the 
notch to correspond to the frequency of the lowly damped mode might lead to instabil- 
ity. For the case under study, the flutter mode remains separated from the notch at 
dynamic pressures near the design point. Furthermore, ground vibration tests indi- 
cate that the frequency of this mode is actually much higher (220 rad/sec) and that a 
substantial amount of structural damping is present. Thus, when the design is 
repeated with the updated structural model, the notch will be widened and occur at a 
frequency much higher than the frequency of flutter. 

Figure 11 presents the Nyquist plot of the system loop transfer function 
obtained by use of the reduced order controller and the evaluation model of the 
plant. It is seen, by comparing with figure 8b, that the reduced order controller 
provides essentially the same stability margins as the full order controller. 
Figure 12 depicts the frequency response of the reduced order controller which is 
essentially the same as that of the full order controller shown in figure 9b. 
Figures 11 and 12 indicate that the order reduction has resulted in little degrada- 
tion in controller performance. 

A high pass filter to remove low frequency signals and additional high 
frequency attenuation will be incorporated in the actual FS implementation. The 

impact of phase changes due to such elements will be considered in future design 
studies. 


Scheduling Law for Improved Off Design Performance 

The closed loop system is unstable at the evaluation points on the desired 
flutter boundary (FS -on) with the controller defined in equation (17). Furthermore, 

F 0 is higher than required at dynamic pressures below that of the design point. 
Performance at off-design points is improved by scheduling the gain as a function of 
dynamic pressure. The schedule is defined in equation (18). 

F = 0.50 F q for q < 400 lb/ft 2 

F - [0.50 + 0.00232(q - 400 )]F q for 400 < q <_ 940 lb/ft 2 (18) 

F = 1.75 F for q > 940 lb/ft 2 
o 

Table 2 shows the performance of the scheduled controller at each of the points indi- 
cated in figure 5. The rms control deflection and rate requirements were computed 
using the frequency domain form (Eq. (2)) of the evaluation model of the plant, the 
scheduled reduced order controller, and a von Karman input gust spectrum having a 12 
fps rms. It was assumed that the measurement noise and fictitious process noise were 
zero. Good phase and gain margins are achieved near the design point. The system is 
stable for cases 7 and 8 which are on the desired flutter boundary (FS-on) of 
figure 5. The control power violations at these two high q cases are not of 

particular concern since flight tests will not be made at these extreme conditions. 
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Optimization of the Robustness of Reduced Order Controllers 

An alternate design approach will now be employed which allows explicit con- 
sideration of design criteria. The approach requires that the form of the control 
law be specified. Consequently, it is particularly applicable in modifying an exist- 
ing control law when small changes occur in the plant or in satisfying design cri- 
teria not fully considered in the previous design. In the discussion to follow, 
robustness design criteria will be emphasized. 

Constrained minimization techniques (Refs. 10-13) can be employed to determine 
the benefits of maximizing the minimum singular value of the return difference trans- 
fer function subject to explicit constraints on gain margins, phase margins, control 
power and controller frequency rolloff characteristics. The nonlinear programming 
algorithm (Refs. 15,16) selected to carry out this optimization requires no gradient 
computations. Some initial results illustrating the procedure follow. 

A reduced order controller of the form of equation (17) was found which maxi- 
mized the minimum singular value of _the return difference transfer function for a 
flight condition at h = 15,000 feet, M = 0.91 (See Fig. 5). Nine of the controller 
coefficients of equation (17) were allowed to vary as indicated in equation (19) 

a^ (s + a 2 ) (s + a_) s z + a,s + a 7 

U FS ~ 2 " * ~~7 T(s) (19) 

s + a 4 s + a 5 s + a 8 s + a g 

where 

692.3 s 2 + 1.031s + 136. 4 2 s 2 + 251 .Is + 314. I 2 , v 

T (s) = _ : (20) 

s + 692.3 s 2 + 2.411s + 136. 9 2 s 2 + 406.1s + 567. 9 2 
was fixed. The following constraints were imposed 

0.25 < a l < 2.25 

40 < a 2 < 120 

150 < a 3 < 300 

50 < a 4 < 175 

2000 < a 5 < 5000 
125 < a 6 < 325 
20,000 < a ? < 35,000 

30 < a g < 90 

6,000 < a g < 12,000 
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> 

> 


N 


c 


(6 a) < 
rms 


($a) < 

rms 


the design 


40° 

6dB 

1 

15° 


740°/sec 

variables were satisfied by making the 


a 


i 



. TT 
sin — 


+ a + a, 

u, J 


J 


( 21 ) 


Here N c is the number of counterclockwise encirclements of the -1 point as 

to varies from 0 to 00 , <(> m is phase margin and F m is gain margin. The 

stability of the closed loop system is determined by evaluating N c * A constrained 
optimization is carried out as follows: 

1. Let c be the vector of inequality constraint functions and define - 


c. = max(0,c,-c t c 0 -c) . 
i 7 i u^ 

2. Define a positive definite diagonal weighting matrix W. 

3. Form the augmented function 

V a '> -IFTT + + 1 o50 |n c -i|- (22) 

min 1 1 

4. Find (a'y)* which minimizes 4>y. 

5. Increase the magnitude of W and repeat steps 3 and 4. 

In the limit as W °>, the solution, if it exists, of (22) will maximize the 
minimum singular value subject to the imposed constraints (Ref. 34). If no solution 
exists, one would have to change the form of the control law or relax the con- 
straints. 

The results obtained in applying this procedure are presented in figure 13 and 
in table 3 for a finite weighting matrix W. The system is stable. The minimum 
singular value is increased 26 percent and gain and phase margin constraints are 
satisfied to within a 2.5 percent tolerance. 
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CONCLUDING REMARKS 


An impleraentable control law has been designed for suppression of flutter in 
the symmetric degrees of freedom for a mathematical model of the DAST ARW-2. Stabil- 
ity margin criteria of ±6dB gain margins and ±45° phase margins are exceeded at the 
design point which is at a dynamic pressure 16.5 percent above the corresponding 
FS-off flutter value. A gain schedule was defined which provided at least ±6dB gain 
and -45 phase margins at off design points examined which had dynamic pressures 
lower than the design point dynamic pressure. The FS law, with gain scheduling, 
resulted in a stable system at points examined on the desired (FS-on) flutter 
boundary having dynamic pressures up to 82 percent above the corresponding open loop 
flutter values. 

Coupling the ISAC program for definition of plant design and evaluation models 
with the ORACLS LQG methodology provided an effective tool for design of full order 
controllers. Addition of process noise at the input allowed stability margin 
criteria to be met. Reduction of the order of the controller from 25th to 9th was 
achieved with minimal sacrifice in controller performance. 

Constrained optimization techniques were successfully applied to maximize the 
robustness characteristics of the reduced order controller at one flight condition. 
The nonlinear programming algorithm used required no gradient computations. Stabil- 
ity the candidate control law at each iteration was determined by computing the 
number of counterclockwise encirclements of the -1 point made by the open loop trans- 
fer function as to varied from 0 to °°. 


New control laws will be designed for both symmetric and antisymmetric degrees 
of freedom as the mathematical model of the ARW-2 is improved. At that time, control 
power and robustness constraints will again be explicitly included in the constrained 
optimization solution for reduced order controllers. 
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le« of Free-Free Symmetric Elastic Modei 



frequency, hz 

i 

7.847 

2 

14.21 

3 

21.72 

4 

30.27 

5 

33.28 

6 

41.10 

7 

47.01 

8 

63.06 

9 

67.22 

in 

78.24 


Table 2.- Performance of scheduled, reduced-order controller. 


case 

q , 

h 

■ 

(6 a)rms 

BfWI 

Cain Margin 

at u 

Phase margin at (■> 

number 

lb/ft 2 

feet 

D 

deg 

deg/ sec 

dB 

rad/sec 

deg 


rad/sec 

1 

547 

18,000 

0.86 



-16.9 

at 

116 

-72.5 

at 

95.6 







13.5 

at 

441 

68.1 

at 

159 

2 

576 

15,000 

0.83 

4.03 

465 

-13.9 

at 

121 

-74.4 

at 

88.2 







12.0 

at 

444 

52.5 

at 

165 

3 

614 

13,400 

0.83 

4.44 

523 

-8.92 

at 

125 

-66.5 

at 

88.0 







10.6 

at 

444 

43.3 

at 

166 

4 

619 

15,000 

0.86 

4.94 

563 

-7.31 

at 

122 

-53.8 

at 

96 







10.8 

at 
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45.2 

at 

162 

5 

670 

13,000 

0.86 

5.41 

629 


Bl 


-47.3 

at 

94.5 
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35.6 

at 

165 

6 

692 

15 ,000 

0.91 

5.76 

599 

-6.08 

at 

122 

-28.2 

at 

86.8 







7.33 

at 

437 

41.8 

at 

174 

7 

912 

8,000 

0.91 

14.2 

1143 

-4.26 

at 

136 

-3.25 

at 

74.7 







2.13 

at 

68.7 

32.9 

at 

203 

8 

938 

4,250 

0.86 

7.90 

837 

-3.43 

at 

150 

-14,2 

at 

75.6 







3.47 

at 

445 

22.9 

at 

190 


Table 3.- Comparison of controller optimized for robustness 
with nominal controller (h-15,000 feet, MKJ.91) 


Function, Constraint , 
and Design Variables 

Optimized Controller 

Nominal Controller 

singular value 

0.587 at 434 rad/sec 

0.465 at 93.4 rad/sec 

(da)rms, deg 

5.23 

5.68 

(ja)rms, deg/sec 

561 

612 

gain margin, dB 

-5.92 at 123 rad/sec 

-4.64 at 122 rad/sec 


7.69 at 433 rad/sec 

8.77 at 437 rad/sec 

phase margin, deg 

-38.8 at 91.2 rad/sec 

-27.1 at 94.6 rad/sec 


39.0 at 177 rad/sec 

40.5 at 164 rad/sec 

. deg/g 

0.543 

0.506 

aj , sec" 1 

106 

98.6 

*2 • «c _1 

234 

213 

■4 , “c -1 

116 

119 

-2 



*5 * “ c 

4529 

4395 

a 6 , see -1 

242 

215 

-2 



*1 , *ec 

26429 

27076 

*8 * meC ~ l 

65.1 

58.5 

-2 


1 

•q. ®ec 

7887 

8915 
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Figure 2. Normal components of six lowest frequency symmetric free-free mode shapes 
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Amplitude , g/deg 



ALTITUDE. FT 



Freq, rad/sec 



Figure 6.- Variation of symmetric characteristic roots with 

altitude for a Mach number of 0.86 (evaluation model). 
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Figure 7.— Closed loop block diagram. 





PHRSE MRRGINS FREQUENCY 
31-33 166.17 

-33.13 gg.91 


GRIN MRRGINS FREQUENCY 
A_. -.676E+01dB 123.91 

B... .GaiE+OldB 257.24 


MINIMUM SINGULAR VALUE 
.463 AT iaa .11 


a) R = 0. 
u 


Imaginary 



PHRSE MARGINS FREQUENCY 
46.25 163.97 

-53.14 95.77 


(SIN MARGINS FREQUENCY 

A. .. -.744E-t01dB 121.68 

B. .. .107E+02d3 440.07 


MINIMUM SINGULAR VALUE 
.690 AT 431.74 


b) R = 0.79 2 deg 2 , 
u 

Figure 8.- Nyquist plots of HG transfer function at the design 
point using full order controller. 

















Imaginary 



PHASE MAR6INS FREQUENCY 
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B— -108E+02dB 441.31 


MINIMUM SINGULAR VALUE 
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Figure II.- Nyquist plot of HG transfer function for reduced 
order controller at the design point. 



Figure 12.- Bode plot of reduced order controller transfer 
function (H) , deg/g, at the design point. 
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Imagincry 


*m m 


-2 


*mF 


Rad 



PHRSE MARGINS FREQUENCY 
40 JO 163.58 

-27.09 94.56 


GAIN MARGINS FREQUENCY 

A. .. -,464E+01dB 122^8 

B. .. ,877E+01dB 436.92 


MINIMUM SINGULAR VALUE 
.463 AT 93.41 


a) nominal, reduced order controller. 


Imaginary 



Rad 


PHASE MARGINS FREQUENCY 
39110 177.45 

-3&80 91.72 


GAIN MARGINS FREQUENCY 
A~. -.592E+01dB 123.48 

B... ,769E+01dB 432.96 


MINIMUM SINGULAR VALUE 

-587 AT 433.78 


b) optimized reduced order controller. 

Figure 13.- Njrquist plot of HG transfer function 
(M = 0.91, h » 15000 ft). 
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